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ABSTRACT 



We consider the Standard Model with an extended scalar sector, and study the 
possible effects of the physics underlying such a model using an effective lagrangian 
parametrization. It is found that certain two photon processes offer windows where 
such heavy interactions might be glimpsed, but the realization of this expectation 
requires enormous experimental precision. 



1. Introduction 



The use of effective lagrangians in parametrizing physics beyond a given the- 
ory has been studied extensively in the recent literature, especially within the 
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context of the Standard Model. The formalism generates a model-independent 
parametrization of the low energy limit of any Green function in terms of a set of 
unknown constants, each of which multiplies a local operator involving only the low 
energy fields; all symmetries are also naturally preserved. Despite the fact that the 
formalism involves, in principle, an infinite number of local operators in the light 
fields, only a finite number of them need to be considered in any given calculation; 
the number of operators which are considered is determined by the required degree 
of accuracy: for higher precision more terms in the effective lagrangian must be 
included, and the number of parameters increases. Any two models of describing 
the heavy physics will generate low energy Green functions whose expressions fit 
the effective lagrangian parametrization (the explicit values of the parameters are, 
of course, model dependent). 

The requirement that the low energy particle content should be the same as in 
the Standard Model is clearly an assumption. One can envisage, for example, the 
case where the Higgs particle is either heavy or absent [4] , or a situation where there 
is an additional relatively light vector boson [5]. In this paper we will consider a 
different modification of the low energy lagrangian for which the low energy scalar 
sector contains two doublets. Such a model, motivated by various supersymmetric 
extensions of the Standard Model [6], will be our low energy theory; the adjective 
"light" will refer to it. 
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Starting from this low energy theory we will investigate the possible effects of 
(heavy) physics underlying it. The purpose of this investigation is dual: on the one 
hand we will determine the difficulties in uncovering heavy physics effects even in 
some very favorable circumstances, which will be of use when and if the minimal 
super symmetric extension of the Standard Model is verified experimentally. On 
the other hand the calculations presented in this paper provide a good illustration 
of the uses of effective operators within the context of renormalization theory. 

Since the light theory is renormalizable and has light scalar excitations (whose 
masses are not protected against 0(A) radiative corrections), the underlying physics 
is expected to be light and decoupling [3]. We will denote by A the scale at which 
the underlying physics becomes apparent, then all observables can be expanded in 
a power series in 1/ A [7] . 

The two doublet model is known to have problems with flavor changing neutral 
currents. We chose to alleviate these difficulties by imposing the usual discrete 
symmetry [8]. We will for simplicity also impose this symmetry on the heavy 
operators (within the framework of this paper this means that the processes which 
violate this symmetry occur at a scale A' ^> A, which is supported by explicit 
computation of the bounds on the scale of process which mediate flavor changing 
neutral currents [9]). We will argue below that the same estimates of the sensitivity 
of a given experiment to A are obtained even when the underlying physics does not 
obey the discrete symmetry. We will also ignore CP violation in the light theory; 

#1 Note that in the process of integrating out the heavy physics some quantities will get 
corrections which do not vanish as A — > oo, but all these contributions can be absorbed in a 
renormalization of the two doublet model parameters and are therefore unobservable. The 
effects of these contributions affect the naturality of the model and will no be considered 
further in this paper. 
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more precisely, we will assume that the scale of CP violation is > A, as discussed 
in section 4. 

A general study of the two doublet model together with all operators of dimen- 
sion > 4 would be both prohibitive and obscure. Fortunately this is not necessary: 
most processes receive light physics contributions (i.e., those generated by the two- 
doublet extension of the Standard Model) which dominate any contribution from 
the heavy physics. Thus we need only concentrate on those processes which are 
suppressed within the two-doublet model and which acquire significant contribu- 
tions due to the inclusion of effective operators. 

As an example of the above situation we will study the two photon decay of 
the CP-odd scalar ao present in this theory, as well as ao production in polarized 
photon colliders. Since this particle has no couplings to the gauge bosons within 
the light theory, the usual contribution to these process are mediated by fermion 
loops. Explicit calculation as well as simple arguments show that these graphs 
are sizable only for the top quark contribution which is proportional to cot (3 
and will be small in the (3 — > n/2 limit. Alternatively we will see that the pro- 
duction of the ao within the light theory in photon collisions (via quark loops) is 
suppressed for a certain combination of photon polarizations. This, however, is not 
the case for the amplitude induced by the effective operators. We can then arrange 
for a suppression of the light-theory production amplitude, with the correspond- 
ing enhancement of the "anomalous" signal, by choosing the photon polarizations 



#2 tan j3 denotes the ratio of the vacuum expectation values, see below. 
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appropriately. 

The organization of the paper is as follows. In the next section we will construct 
the relevant operators for the above processes. The calculation of the corresponding 
amplitudes together with some illustrative numerical results is presented in section 
3. Section 4 contains some parting comments. Details of the calculations are 
summarized in the appendix. 



2. The Lagrangian 

We begin with a brief description of the model starting with the renormalizable 
interactions; later we turn to the study of heavy physics effects using effective 
operators. 

The scalar sector contains two doublets 4> a (a = 1,2), while the fermionic and 
gauge particle contents are identical to those in the Standard Model. As mentioned 
in the introduction we will impose a Z2 symmetry on the low energy sector in order 
to suppress flavor changing neutral currents. Under this symmetry <pi and all right 
handed fermions are odd, while the remaining particles are even. We will use the 
conventions of Refs. 9 and 10. 

The light scalar potential (that is, the terms in the potential of dimension < 4) 

#3 The constants multiplying the effective operators are estimated using standard arguments 
(see below); it is of course possible for some of these constants to be suppressed due to 
unknown effects, we will comment on this possibility later on. 
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is given by 



V =Ai - v^ 2 + A 2 (0 2 02 - *|) 2 

+ A 3 (0l01 + 0^02 - V 2 ) + A 4 (0l0l) (0^02 



^102 



(2.1) 



+ A 5 



Re 



^102 



V\V 2 



+ A 6 



Im 



^102 



In terms of the mass eigenstates the scalar fields are 



1 i 
i>i = vi + -j=(c a H - s a h ) + -^(cpGo - spao) 

1 % 
L >2 = v 2 + -^{s a H + c a h ) + -j=(spGQ + cpao) 



(2.2) 



0+ = c p G+ - s p H+ 
s p G + + c p H+ 
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where 



0a 



<f>t 



a = 1,2 



(2.3) 



We used the usual definitions tan/5 = vifvx, v z 



cos a, etc. 



The fields Go, G + and G correspond to the would-be Goldstone bosons (with 
G~ = (G + )^), while the other scalar fields are physical. 

In the calculations below we will adopt the following notation and conventions 
(see Ref. 9): £ and q denote the left-handed lepton and quark doublet respectively, 
e, u and d denote the corresponding right-handed SU{2) singlets, (of particular 
interest to us will be the case where u denotes the right-handed top quark and 
q the left-handed top-bottom doublet). The SU{2) gauge fields are denoted by 
(where / is a weak isospin index) and the U(l) gauge field is labelled B^, 



the corresponding field strengths are W^ v and i?^ ; finally r denote the Pauli 
matrices and (fik = iT 2 (fi* k . The covariant derivatives are = d^—^^W^—^-YB^ 
where Y is the hypercharge and g, g' the gauge couplings corresponding to SU(2) L 
and U(l)y respectively. 

With these preliminaries we can now discuss the possible effects of heavy 
physics on low energy observables (where "low energy" refers to scales below that 
of symmetry breaking in the above two-doublet model), concentrating on the cor- 
responding modifications to the coupling of the CP odd scalar ao to two photons. 
We will assume that the underlying physics is decoupling and weakly coupled, 
and denote by A its characteristic scale (the scale at which it becomes directly 
observable). We only consider the weakly coupled case since for strongly coupled 
underlying physics it is difficult to maintain the scalars <p a naturally light (com- 
pared to A) without considerably modifying the low energy spectrum [3]. The 
decoupling theorem [7] then implies that all the effects of heavy excitations at low 
energies can be parametrized in a model independent manner by a series of local 
operators which satisfy the same symmetries as the low energy theory, and whose 
coefficients are suppressed by the appropriate power of A. It is of course possible 
for the underlying theory to have several mass scales, if this is the case A will 
denote the smallest of these scales.^ 5 

The construction of the effective lagrangian closely parallels the one described 
in Ref.9. The terms of dimension four are the ones in the two-doublet model; 
there are no terms of dimension five due to Lorentz and gauge symmetries (right- 

#4 is the full non-Abelian expression for the curvature. 

#5 This implies that a suppression factor <~ (A/A') 2 may be present in some of the operators. 
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handed neutrinos are assumed to be absent). Therefore all the effects of the heavy 
excitations will be suppressed by at least two powers of A. The simplest way 
of obtaining these operators is to take the list provided in Ref.9 and replace the 

#6 

Standard Model doublet by either <\>\ or 02 

Each of the effective operators appears in the effective lagrangian multiplied 
by an undetermined coefficient. This apparently implies that there is no way 
of extracting from this formalism quantitative estimates. Fortunately this is not 
the case: having taken the underlying physics to be weakly coupled, the order of 
magnitude of the coefficient of an operator is determined by the type of graph in 
the underlying theory which generates it. If an operator is generated at tree level 
by the heavy dynamics, it will appear with a coefficient ~ h/A 2 where h denotes 
a product of coupling constants in the underlying theory; these by assumption are 
< 1, so that we also expect h<l. Loop-generated operators will acquire an extra 
suppression factor, so that their coefficient are of the form h'/^nA) 2 with h' <1. 
The additional loop factor ~ \/(Aix) 2 insures that the loop-generated operators 
are subdominant. We will assume that the heavy physics is described by a gauge 
theory and, for quantitative estimates will take the corresponding gauge coupling 
to be ~ 1; one should keep in mind that the results thus obtained represent an 
upper bound and could be suppressed by a significant amount. 

The standard (two-doublet) contribution to the ao-7-7 coupling occurs at one 
loop [10]. The only dimension six operators which generate tree-level contributions 

#6 We have verified that this procedure generates all the operators relevant for the calculations 
below, despite the dfact that the equations of motion were used in Ref. 9 to eliminate some 
operators. 



8 



to this vertex are 

Ctreel = {B^f Im (<j>\<fo) , O t reel = B^B^ Im (0{02j , 

2 (2-4) 
Otree2 = Im (^ 2 ) , a trcc2 = Im (0l02) ; 

(together with their hermitian-conjugate counterparts), which are loop-generated 
[11]: there are no dimension six, tree-level-generated operators that contribute to 
the ao-7-7 vertex at tree level. Note that these operators violate the Z2 symmetry. 
A consistent loop expansion then requires that we include all tree-level contribu- 
tions containing one (2.4) insertion. We must also include all one-loop graphs 
containing one O insertion, provided this operator is generated at tree level. The 
one loop graphs containing insertions of loop-generated operators can be ignored; 
this considerably reduces the number of terms that need to be considered. 

At this point we can either assume that the underlying dynamics respects the 
discrete symmetry or not. The precise results differ depending on which of these 
cases is realized; note however that in either case all contributions to the vertex of 
interest occur at one loop; either through tree-level-generated operators in one loop 
graphs, or through loop-generated operators in tree level graphs. In this paper we 
are only interested in estimating the effects of the underlying physics and for this it 
is sufficient, in view of the previous comments, to consider the situation where the 
underlying physics obeys the discrete symmetry. This assumption has the added 
advantage of simplifying the calculations. 

When the physics underlying the two doublet model is assumed to respect the 
discrete symmetry there are no tree-level contributions to the ao-7-7 vertex since 
the operators in (2.4) are forbidden. It follows that the dimension six modifications 
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of the low energy lagrangian contribute to the processes under consideration only 
through (explicit) one loop graphs; as mentioned above only tree-level-generated 
operators need be considered. 

Due to the absence of tree level interactions, the loop contributions must be 
finite. Indeed, any divergence can be associated with a local operator respecting 
the symmetries of the model; the effective lagrangian already contains all such 
operators, and all divergences can be absorbed in a renormalization of the effec- 
tive lagrangian coefficients [12]. Since there is no tree level operator (when the Z2 
symmetry is imposed) containing an ao-7-7 term, it follows that the correspond- 
ing three point function must be finite, else it would be impossible to absorb its 
divergence. The cancellation of divergences, together with electromagnetic gauge 
invariance, provide important checks on the algebra. Imposing the discrete sym- 
metry on the effective operators also implies that there is no dependence (to one 
loop) on the renormalization scale. 

The effective lagrangian takes the form 



where the summation over % runs over the set of (tree-level-generated) operators 
given in the appendix. The unknown constants y% depend on the underlying physics 
and, based on the discussion above, are expected to be ~ 1. 

The most important terms in (2.5) are generated by the operators 




(2.5) 



O uip -211 = 



(4<h) (Wi) , 0% = [d, ; (2.6) 



which dominate the heavy physics contributions to the ao-7-7 vertex (in i, j, k, I 
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1, 2 and only those terms with an even number of <\>\ factors are considered due to 
the discrete symmetry). The corresponding contributions are not suppressed in the 
limit of large tan/5; nor are they suppressed by (explicit) small fermions masses. . 



3. Calculation of the amplitudes 

Having described the lagrangian we can use it to evaluate the three point 
function for two photons and one ao The photons will be labelled by the subscripts 
1 and 2, the corresponding momenta and polarization vectors will be denoted by 
k a and e a (a = 1,2) respectively. We will assume that all (external) particles are 
on shell. 

Due to electromagnetic gauge invariance there are two possible tensorial struc- 
tures allowed for the amplitude: 

T>i = (ki ■ k 2 )(ei ■ e 2 ) - (h ■ e 2 )(/c2 • ei), V2 = «e a /3^ e i e 2 k i ^2 ( 3 -l) 
correspondingly the amplitude can be written as 

Aao-Tl = M1V1 + M 2 V 2 . (3.2) 

The light-physics amplitude generated via quark loops contains only Vi- In 
contrast the O induced terms can be proportional to either quantity, in particular 
all bosonic loops give results proportional to V\. This implies that the decay ao — * 

^7 It is of course possible for the underlying dynamics to generate coefficients which are small 
for large tan/3 and/or contain small Yukawa couplings, we will not consider this possibility 
here. 
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77 w iH have a modified width and angular distribution due to the presence of the 
operators O. Unfortunately the branching ratio is so small that the modification 
of the angular distribution is unobservable (at least in the forseeable future). 

In contrast, the presence of terms proportional to V\ in the amplitude can 
have important effects in ao production in photon colliders. These machines will 
(probably) allow the possibility of polarizing the photons whence the terms oc T>2 
can be suppressed by the appropriate choice of initial polarizations thus enhancing 
the heavy physics contributions. 

The calculation of the amplitude involves both one particle irreducible and one 
particle reducible diagrams. The latter arise due to the ao — ho and aQ — Ho mixings 
induced by operators such as the second one in (2.6). The main contributions to 
the 1PI graphs come from the first operator in (2.6), which also becomes the 
overall dominant term when a is close to zero or to n/2 (provided tan/5 is large). 
For intermediate values of a the contributions from both operators in (2.6) are 
comparable. 7 ^ 8 

The total amplitude derived from (2.5) is 



where Ai denotes the contribution to the ao-7-7 (on-shell) three point function 
generated by Oi (without the factor m/A 2 ). From Ai ^77 the two-photon decay 
width of the ao or its production cross section in photon collisions can be evaluated 

#8 The 1PR graphs will in general exhibit resonances at the mass of the ho and Ho. The 
previous comments correspond to energies outside these regions. 




(3.3) 
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directly. Note that the terms containing V\ and Vi will not interfere in polarization- 
averaged widths or cross-sections. The quantities Ai are displayed in the appendix. 
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T(a -» 77) X 10 5 ( GeV) r(a -» 77) X 10 6 ( GeV) 




™a (GeV) 



(c) (d) 



Figure 1. Width and branching ratios for the decay ao — ► 77 as a function of m O0 for two 
values of /?. The solid curves correspond to the theory with effective operators for 
A = 1 TeV, the dashed curves give the pure two-doublet model results. We chose 
m hg = 450 GeV, m Ho = 550 GeV, m H+ = 500 GeV; m top = 174 GeV, a = 45°. 
The resonance peaks are produced by the 1PR diagrams. 



For generic values of the parameters in the standard two-doublet lagrangian 
the contributions from the Oi are small due to the presence of the 1/A 2 factor. This 
suppression can be countered in two ways. For ao production we can use the photon 
polarization to suppress the light contribution and bring the anomalous effects to 
the foreground. For the two-photon decay of the ao the light contributions are 
proportional to either mbtan/3 or m t cot (3 and are suppressed when \fl — ir/2\ is 
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similar to zero (though larger than ~ m^/m^). The only terms which can compete 
with the light contributions are those which do not vanish as (3 — > n/2 and are 
not suppressed by a small mass factor. These constraints are satisfied only by the 
operator (2.6). 

Using the results of the appendix we evaluated the ao — > 77 width and branch- 
ing ratio for several illustrative choices of the masses and couplings, the results are 

#9 

presented in Fig. 1 

r(a -» 77)/r< SM >(# -» 77) r( ao -» 77)/r< SM »(^ -» 77) 




400 600 S00 1000 



m O0 (GeV) 

(c) (d) 

Figure 2. (a),(b): Ratio of the total ao — ► 77 decay width to the decay width of a 
Standard Model Higgs particle of the same mass into two photons. (c),(d) Total 
number of yearly events at the LHC (CM. energy^ 7TeV, luminosity 100/fb). The 
solid curves correspond to the theory with effective operators for A = 1 TeV, the 
dashed curves give the pure two-doublet model results. We chose m^ = 450 GeV, 
m Ho = 550 GeV, m H+ = 500 GeV and A = 1 TeV; m top = 174 GeV, a = 45°. 



#9 These results are presented for the choice of constants j/j specified in the Appendix. 
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As can be seen from this figure there is a ~ 50% increase from the light-physics 
values for this decay width; despite this one has to face the complications produced 
by the smallness of the branching ratio, being, in the most favorable case ~ 1CP 7 
and about 3% of the two-photon width of a Standard Model Higgs boson of the 
same mass, as shown in Fig. 2a and 2b. The number of events, showed in 2c, d is 
marginal for one year of LHC running, but, mainly due to the resonance enhance- 
ment, deviations from the two doublet model would be observable for moderate 
values of m ao . 

This implies that when (and if) the ao is discovered at the LHC, its branching 
ratio into two photons could prove a sensitive reaction in which to study physics 
beyond this excitation but the experimental sensitivity required is enormous. We 
have verified that all heavy physics effects are completely obscured when tan (3 is 
smaller than ~ 5; the results are also insensitive to the mass of the charged Higgs. 

We also investigated the behaviour of the two-photon branching ratio as a 
function of A; the results are presented in Fig. 3. Based on this graph we conclude 
that this reaction will be able to probe physics up to ~ 1 TeV. 

B(ao -> 77) x 10 s 




400 600 800 1000 

m ao (GeV) 



Figure 3. Branching ratio for the decay ao-7-7 for A = 1 TeV (solid curve), A = 2 TeV 
(dashed curve) and A = 00 (dotted curve). We chose mh = 450 GeV, mjy = 
550 GeV, m H + = 500 GeV and tan/3 = 30; m top = 174 GeV, a = 45°. 
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Now we consider ao production in photon-photon collisions. From general 
considerations the complete amplitude is of the form *4 77 _,. ao = M.{P\ + 
We will assume that both initial photons have the same polarization matrix. In 
this case we obtain 



\A I 2 



£ [iMii 2 (i+^-e + ) + i^ 2 | 2 (i-e+) 



(3.4) 



where £ ± = (£1, ±£2, £3) and £j are the Stokes parameters for the photons (assumed 
the same for both). 

The ratio (1 — £+)/(l — £ is independent of the polarization of the intial 

photons and depends only on the polarization of the initial electrons (for perfectly 
polarized electrons £ + = 1). As an illustration we will assume that the initial 
electrons are 70% longitudinally polarized, the electron and photon energies are 
taken equal to 500 GeV and 1.17 eV. In this case [13] (1 -£+ )/(l -£ + ) ~ 0.86 
(and becomes ~ 0.7 at 90% polarization) 



tan/3 = 10 




00 1000 
m ao (GeV) 



(a) 



Figure 4. Total ao production cross section averaged over an energy bin equal to the total 
width of the ao- The solid curves correspond to the theory with effective operators 
for A = 1 TeV, the dashed curves give the pure two-doublet model results. We chose 
m ho = 450 GeV, m Ho = 550 GeV, m H+ = 500 GeV; m top = 174 GcV, a = 45°. 
The electron polarization is 70% . 



(7(77 -> ao) (fb) 




400 



600 



(b) 



#10 For a more significant reduction exquisite degrees of polarization are required: 99% for 
(!-£)/(! 0.25 
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These results imply that, even though it is theoretically possible to eliminate 
the light-physics contribution to this cross section for any value of f3, experimental 
constraints on the initial electron's degree of polarization prevent the realization 
of this possibility. The expression for the cross section is simply (in the center of 
mass system) 

dcr=— — — \A^^ ao \ 2 S (y/s - m ao ) . (3.5) 

ZS 7Jl ao 

In figure 4 we plot the value of this expression averaged over the width of the 
ao, namely 



7T | .A-y-y—^ao | 

z " L a 1 do 



(3.6) 



To illustrate the significance of this result we consider the ratio of a to the same 
quantity when A — > oo. As can be seen from Fig. 5, the deviations are important 
for A = 1 TeV. 



<t(77 -> a )/<T 2 doublet(77 -* a o) 




400 600 800 1000 

m ao ( GeV) 



Figure 5. Ratio of the total to light cross sections. The solid curve correspond to tan/3 = 
30, the dashed curve to tan/3 = 10. We chose m/, = 450 GeV, iriH — 550 GeV, 
m H + = 500 GeV and A = 1 TeV; m t0 p = 174 GeV, a = 45°. The electron 
polarization is 70% . 



The number of events to be generated at a 77 collider with 10/fb/yr luminosity 
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is presented in Fig. 6. As in the previous cases, deviations from the usual two- 
doublet model are noticeable for large values of tan (5 only. 



tan/3 = 10 




00 1000 
m O0 (GeV) 



(a) 




400 



600 



GO 



Figure 6. Number of a events produced in 77 colissions for two values of tan/3. We 
chose m ho = 450 GeV, m Ho — 550 GeV, m ao = m H + = 500 GeV; m top = 174 GeV, 
a = 45°. The electron polarization is 70% and the luminosity equals 10/fb. 



AAf/VN AAf/VTJ 




1500 2000 2500 3000 

A (GeV) 



(a) (b) 

Figure 7. Statistical significance of the deviations from the light cross section, solid curve 
tan/3 = 30, dashed curve tan/3 = 10. We chose nrih = 450 GeV, to# = 550 GeV, 
m ao = m H + = 500 GeV; m top = 174 GeV, a = 45°. The electron polarization is 
70% and the luminosity equals 10/fb. 

Finally, in order to estimate the statistical significance of the deviations from 
the light-physics result, we evaluate A/", the total number of events and AM, the 
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total number of events minus the number of events present when A — > oo. Then 
AJ\f I a/AT is a measure of the statistical significance of the deviations from the 
light-physics results and is plotted in Fig. 7 for a 10/fb luminosity. As can be seen 
from this plot there are significant deviations from the light-physics results up to 
a few TeV. 

4. Conclusions 

We have employed the method of effective lagrangians in evaluating the possible 
deviations from the two-doublet extension of the standard model for reactions 
containing the CP odd excitation ao and two photons. The calculations illustrate 
the fact that effective lagrangians can be used in loop calculations. 

The magnitude of the resets can be estimated reliably within a set of general 
scenarios. For example, we must examine whether CP violation operators are gen- 
erated at the same scale as the CP conserving ones, or whether the underlying 
physics is weakly or strongly coupled. One of the advantages of the effective la- 
grangian approach is the ease with which these various possibilities are identified 
and studied. 

The cases which we studied in detail were the two-photon processes related to 
the CP-odd scalar ao- We found important deviations stemming from the anoma- 
lous couplings when the scale of new physics is moderate large (a few TeV) provided 
tan/?> 20. There are several problems with the type of processes considered stem- 
ming mainly from the small branching ratios and cross sections. For the LHC, 
the event number is so reduced (see Fig. 2) that the observation of the effects 
generated by the physics beyond the two-doublet model is very unlikely (except 
for a 10-year study). In constrast, for a 77 collider with a luminosity of ~10/fb, 
a careful study will uncover some effects generated by the physics beyond the two 
doublet model, or at least place significant bounds on the scale of this type of new 
physics (see Figs. 6,7.). 
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It was mentioned in the introduction (as evident from the form of the potential) 
that no CP violation effects were included in the light theory. In contrast, several 
effective operators (such as 0^' kl ) violate CP. This situation is not inconsistent. 
In fact, CP violations in the potential can be pictured as ho — ao and Ho — clq 
mixings of the same type as those induced by 0^ 1,3 \ Therefore the presence of CP 
violating terms in the potential will not affect the conclusions as long as their order 
of magnitude is the same as that of the mixing terms generated by the effective 
operators. More precisely, we have assumed that CP violation is generated at a 
scale > A. 

The numerical results were obtained not by the most optimistic choice of pa- 
rameters, but by mimicking the possible presence of cancellations among various 
graphs. As noted in the appendix there is still the possibility of having further 
suppression factors, in which case the new physics effects will be too small to be 
observed. On the other hand there could be some enhancements, in which case our 
results would constitute a lower bound on the new physics effects. 

The results presented in this paper are to be compared to the ones derived in 
the minimal Standard Model where the (CP-even) #77 interaction is generated 
through one-loop effects of charged fermions and W gauge bosons. The contri- 
butions from the W boson and the top quark are dominant, although the latter 
is only marginally important and tends to cancel the first one partially. When 
the contribution of dimension six operators is considered, the virtual effects may 
enhance the Standard Model width of this rare decay mode up to one order of mag- 
nitude. 14 ' 15 Recently it was pointed out that tree-level generated bosonic operators 
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of dimension eight may also enhance the Standard Model result. 

The conclusions of this paper are, by necessity, speculative. The ao has not 
been found as of yet and so the study of its decays lies in the future. Still, based 
on the above calculations it is clear that if this excitation is observed, the detailed 
study of it's interactions with two photons may open a window into new physics. 
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APPENDIX 

In this appendix we present the full list of operators contributing to the ao-7-7 
three point function and the corresponding amplitudes. 

The list of (tree-level-generated) operators which contribute at one loop is the 
following (i, k, l,m,n — 1, 2) 



0&p;ijk = (4>t<i>j)(te<j>k) ; 
Odip-ijk = (4>t<i>j)(qd(j>k) I 

M = i(*tD ll (l> k )(f>ff); f = e,u,d; 
0f kF = * (4 D ^k) {F-fF) ; F = £,q; 



(3) 



0<p;ijklmn 
;ijkl ~ 



1 I o k . ;' !),,< 



■ 95 



4h) (4b) {<tL<t>n) ; 



o 



(1) 

ijkl 



d, 



(3) _ 

w ijkl - 



4D^ 3 ] \{D^ k )Ul 



(A.l) 



where only operators consistent with the discrete symmetry are retained, for exam- 
ple O eip] ijk is considered only with i+j + k odd, O^ijkimn when i+j + k + l + m+n 
even, etc. As mentioned above, we have adopted the notation of Ref. 9. It is worth 
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pointing out that the operator 0\- kl is slightly different from the one used in Eq. 
(3.14) of Ref.9, the difference between this expression and ours is easily seen to 
be, using the equations of motion, a linear combination of 0(p-ijklm 

and Oq^ 

;ijkli 

therefore both expressions are equivalent [16]. The one we chose enormously sim- 
plifies the calculation of the corresponding amplitude (reducing contributions from 
19 to 3 graphs). 

(1 3) 

The operators O^jkimn, Od(p;ijkl an d 0^ kl , generate quadratic terms which 
mix ao with the scalars, h° and H° (this is a consequence of the fact that these 
operators violate CP). Therefore these operators produce two kind of graphs in the 
unitary gauge: 1PI diagrams with H + in the loop, and 1PR diagrams where the 
ao converts into a h° or H° via the above mixing terms, and then the hP or H° 
decay into two photons via H + , W or fermion loops. The amplitude corresponding 
to the above O via 1PR graphs will be denotes by a "(lPi?)" superscript (no 1PI 
superscript will be used in the other amplitudes in order to simplify the notation) 

We denote by Ai the value of the on-shell three point function generated by 
Oi (without the factor m/A 2 ), straightforward evaluation of the relevant Feynman 
rules and graphs yields, for the 1PI contributions containing bosons in the loop, 
generated by O vijk i m ,O dvijk i and { ^ w (a = 1,3) yield (/ stands for the right- 
handed fermions f = e,u, d) 




(A.2) 
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e m w s w 




2— |i/_i(m jff+ ) + 1 Vi; 

. m a J 



a 



= 1,3 



4tt 2 



where mjj+ denotes the mass and 



l 




n > -I. 



(A.3) 
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The coefficients in (A. 2) are 



1.(1) -(_\k 2. l(1) _ (_)k r 2. l(1) - -Ir-o 

"llfcfc — ^ J 6 /3' "22fcfc — ^ J C /3' "1212,1221,2112,2121 _ 2 C 2/3- 

(A.4) 

"2211, 2121 — + "1122,1212 — x ' "1111,2222,1221,2112 — u - 



The 1PI contributions with fermions in the loop are 



_ spb{jk m l s lQ 2 f m f \ 



2 v / 2 7T 2 m 2 
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2pQ}(-) 



k+1 



8tt 2 



0,() 



m. 
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0() 



^(1,3) _ es w rawS2/3Q/( ) fc+1 y-y 



8tt 2 
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2 ^f / -iK) + 1 7,2 
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(A.5) 



where 



ie,d 2 2 2 -i , 2 t m _ 2 2 i 2 

°111,221,212,122 - ~ C /3> _S /3' _S /3' 1 + C /3' °112,121,211 — c /3 ? c /3> -i _ s /3 



4»» 



-1 for a = 3, / = u (A.6) 
+1 otherwise 



Finally the 1PR contributions are generated, as mentioned above, by the con- 
tact terms in some operators which induce ao — h° and ao — H° mixings, the h° 
and H° then decay into two photons via fermion, W or H + loops (in the unitary 
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gauge). The results can be extracted from [17, 10], the amplitude is (0 = h°, H°) 



m{4>) es w m^ 

m a - m l 2n2 
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1 + 



2m 
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I-i(m H +) 
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where Nf denotes the number of colors for fermion /. The no n- vanishing constants 
m{4>) are 
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1212,2121 
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/l0\(3) _ /l0\(3) 1 
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/tr0\(3) _ /rrCh(3) 1 
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^(^°)d</>;1212 = -V(h°)d<p;2121 = 2 C a-/3 _ ^ c 2/?C a _/ 



?7(# )<9</?;1212 = -^(-^°)av3;2121 = ^a - / 8 ~~ | c 2/3Sa-/3 
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f h 0^ (h 0s S l m l i l ~ (-) fcc 2/3) 

r}{ri ) v -kk\2Y2 = ~r]{h ) v -kk2\2\ = 7—— 9 



47ram2 



V(H°)^-kkl212 — - r l(H°) Lp; kk2121 
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C a /C[3 



The total contribution is then 



•A — -^.standard + ^ 2 ( ^ ] ^ ] Z/i^-i ^ J 

V i i J 



(A.10) 



where the constants yi contain all couplings from the underlying theory which 
generates the operators 0{. From this expression the decay width of the ao or its 
production cross section can be evaluated directly. Note that the terms containing 
V\ and V2 will not interfere in polarization averaged widths or cross-sections. 

It must be kept in mind that all the coefficients will contain, in general, some 
coupling constants of the underlying theory which are small ( ^ 1) by our assump- 
tion of it being a weakly coupled theory. For illustrative purposes we will choose 
\yi\ = 1 and we will neglect all anomalous contributions which are proportional 
to a fermion mass other than the top. In order to simulate possible cancellations 
between the various contributions of the same type (produces by constants m of 
opposite signs) we will use an averaging procedure such as 



^ ^ ytip;ijkAt(p-ijk 
y/s=112,121,211 



ijk=112, 121,211 



(A.ll) 



This simple procedure implies several cancellations, for example, the contributions 



2G 



from Afaf and A^f cancel. 

We emphasize that this example is presented for illustrative purposes only, 
nonetheless we expect the results to provide good semi-quantitative estimates. 
Should there be fewer cancellations the anomalous signal would be enhanced; in 
this case the results presented constitute a lower bound. On the other hand it is 
possible for the constants yi could be suppressed by unknown effects in which case 
our results would over-estimate the anomalous effects. 
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FIGURE CAPTIONS 



1) Width and branching ratios for the decay ao — > 77 as a function of m ao for 
two values of (3. The solid curves correspond to the theory with effective 
operators for A = 1 TeV, the dashed curves give the pure two-doublet model 
results. We chose = 450 GeV, m# = 550 GeV, m#+ = 500 GeV; 
""itop = 174 GeV, a = 45°. The resonance peaks are produced by the 1PR 
diagrams. 

2) (a),(b): Ratio of the total ao — > 77 decay width to the decay width of a 
Standard Model Higgs particle of the same mass into two photons. (c),(d) 
Total number of yearly events at the LHC (CM. energy= 7TeV, luminosity 
100/fb). The solid curves correspond to the theory with effective operators 
for A = 1 TeV, the dashed curves give the pure two-doublet model results. 
We chose m ho = 450 GeV, m Ho = 550 GeV, m H+ = 500 GeV and A = 1 TeV; 
mtop = 174 GeV, a = 45°. 

3) Branching ratio for the decay ao-7-7 for A = 1 TeV (solid curve), A = 2 TeV 
(dashed curve) and A = 00 (dotted curve). We chose = 450 GeV, 
m Ho = 550 GeV, m H + = 500 GeV and tan/5 = 30; m top = 17 '4 GeV, 
a = 45°. 

4) Total ao production cross section averaged over an energy bin equal to the 
total width of the ao- The solid curves correspond to the theory with effective 
operators for A = 1 TeV, the dashed curves give the pure two-doublet model 
results. We chose = 450 GeV, mn = 550 GeV, m#+ = 500 GeV; 
m>top = 174 GeV, a = 45°. The electron polarization is 70%. 

5) Ratio of the total to light cross sections. The solid curve correspond to 
tan/3 = 30, the dashed curve to tan/3 = 10. We chose = 450 GeV, 
m Ho = 550 GeV, m H+ = 500 GeV and A = 1 TeV; m top = 174 GeV, 
a = 45°. The electron polarization is 70% . 

6) Number of ao events produced in 77 colissions for two values of tan/3. We 
chose m/j = 450 GeV, m# = 550 GeV, m ao = m H + = 500 GeV; m top = 
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174 GeV, a = 45°. The electron polarization is 70% and the luminosity 
equals 10/fb. 

7) Statistical significance of the deviations from the light cross section, solid 
curve tan/3 = 30, dashed curve tan/3 = 10. We chose rrih = 450 GeV, 
m Ho = 550 GeV, m ao = m H+ = 500 GeV; m t0 p = 174 GeV, a = 45°. The 
electron polarization is 70% and the luminosity equals 10/fb. 
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